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l — ij In this paper the Martin-Siggia-Rose (MSR) functional integral represen- 

tation is used for the study of the Langevin dynamics of a polymer melt 

m : 

t— I , in terms of collective variables: mass density and response field density. 

o ■ 

t^. ' The resulting generating functional (GF) takes into account fluctuations 

as : 

around the random phase approximation (RPA) up to an arbitrary order. 
The set of equations for the correlation and response functions is derived. 

: 

It is generally shown that for cases whenever the fluctuation-dissipation 

> '■ 

• rH , theorem (FDT) holds we arrive at equations similar to those derived by 

><■ 

H 

Mori-Zwanzig. The case when FDT in the glassy phase is violated is also 
qualitatively considered and it is shown that this results in a smearing out 
of the ideal glass transition. The memory kernel is specified for the ideal 
glass transition as a sum of all water-melon diagrams. For the Gaussian 
chain model the explicit expression for the memory kernel was obtained 
and discussed in a qualitative link to the mode-coupling equation. 



I. INTRODUCTION 



Over the last decade the glass transition theory as well as corresponding experiments 
were strongly influenced by the mode coupling approach (MCA) This approach 

based on the Mori-Zwanzig projection formalism || specified for the two slow variables: 
mass density and longitudinal momentum density. The subsequent projection of the 
random forces, which are involved in the memory kernel, onto products of the two densi- 
ties and the factorization of the resulting 4-point correlators yields the closed non-linear 
equation for the density time correlation function 0(k, t). The bifurcation analysis of 
this equation |TJ shows that at some critical values of the coupling constants and con- 
trol parameters (like temperature T, density n, etc.) the non- vanishing long time limit 
0(k, t — > oo) = /(k) arises. This indicates the occurance of a non ergodic (glass) state. 

A complementary approach based on the non-linear fluctuating hydrodynamic (NFH) 
was developed in [[§-[{[]. The authors start here from a set of stochastic equations for the 
mass density and the momentum density. Then by using the renormalized perturbation 
theory and one-loop approximation we derive basically the same equations as in the MCA. 

As distinct from the perturbative treatment the numerical solution of these equations 
is (with an accuracy of numerical errors) exact. Therefore, a comparison of the numer- 
ical solution with the perturbative results provides an estimation of the validity of the 
approximations made in the analytic studies. The results of these calculations show 
that in spite of the fact that the time-dependent density correlation functions are slightly 
stretched the very important features, such as two-step relaxation regime, can not be 
obtained. This regime was reproduced in the numerical study of the Langevin equations 
of the system with a free-energy functional of the Ramakrishnan-Yussouf (RY) form ||10|| . 
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The RY free energy functional provides a large number of glassy local minima, but from 
the results of ref. it is still not clear whether the observed two-step relaxation regime 



arises from non-linearities of density fluctuations in the liquid or from transitions between 
different glassy minima. Moreover, it is not evident how reliable the results of the NFH 
are near the main peak of the static structure factor S(k), i.e. k rs k . In this area the 
density fluctuations are rather strong and the correlator of the noise fields in ref. M,TU 



is 



taken artificially suppressed (the dimensionless parameter A = 10~ 4 ) otherwise the density 
could be even negative. 

In this situation it is very instructive to start from the opposite, in some sense, to the 
NFH limit. In contrast to the NFH the random phase approximation (RPA) describes 
the modification of the behavior of free-particles by effective interactions. The dynamical 
version of the RPA was used mainly for the description of the dynamic spectrum of simple 



liquids [11] as well as for polymer melts and solutions [12-15 . As have been shown in 



ref. [|TT] the basic defect of the RPA in the context of the glass transition problem is the 
absence of a central peak which is determined by the corresponding slow dynamic IjJH]. 
The reason of this is evident: the RPA completely ignores intermolecular collisions which 
only bring the system into a state of local thermodynamic equilibrium and eventually 
assure the hydrodynamic regime. 

In the RPA the particles are free or only weakly interacting. Since these effects in- 
variably dominate at sufficiently short wavelengths, such representation might be suitable 
for description at large wavelengths, k m ko. To assure the feedback mechanism however, 
which is responsible for the glass transition or the microphase separation in block 



copolymers fi6HT9f, approaches beyond RPA are necessary. 



For the static case it is such extension which was carried out in ref. [18 JOS. In ref 



20fl , by using a nonperturbative Hartree approximation, we have been able to derive a 
generalized Rouse equation for a tagged chain in a melt. The freezing process of the Rouse 
modes of the test chain was sequentially considered. In the present paper we emphasizes 
on a systematic way of taking into account density fluctuations in a homopolymer melt 
with respect to the dynamic RPA. In doing so the glass transition dynamics will be 
of our prime interest. This appears as a fundamental problem in polymer physics. The 
RPA is well known in describing several collective phenomena in interacting polymer 
systems quite well. The phase behavior (statics and dynamics) of polymer mixtures, 
block copolymer melts can be understood very well. The theoretical description of freezing 
processes, however, are certainly beyond the random phase approximations. Moreover, 
such freezing processes cannot be of perturbative nature. Interactions become strong 
and dominant on short length scales. Thus we must use methods that go systematically 
beyond the classical RPA in polymer physics. In doing so, the glass transition dynamics 
in the polymer melts will be of our prime interest. We must mention one important point 
here. So far we restrict ourselves to the cases of low molecular weight melts. This is to 
avoid additional complications with reptation dynamics for melts consisting of chains with 
a large degree of polymerization. Our point in this paper is thus to develop a method 
which allows to study the freezing of an ensemble of "Rouse chains" with a degree of 
polymerization below the critical molecular weight N c . 

The paper is organized as follows. In section II the Langevin dynamics of a homopoly- 
mer melt is treated by using the Martin-Siggia-Rose (MSR) generating functional (GF) 



method fl21| , |22|| . The effective action is represented in terms of the two collective variables, 
the mass density and the response field density. In section III the equations of motion 
for the time- dependent correlation function and the response function are derived. It is 



shown that in the regime when the fluctuation dissipation theorem (FDT) is valid these 
two equations are reduced to one, having the form of a MCA-equation, which yields under 
certain conditions an ideal glass transition. The case when the FDT is violated is also 
briefly considered. It is shown that by replacing the usual FDT by the assumption of a 
Quasi-FDT (QFDT) this leads to a smearing out of the ideal glass transition. In section 
IV the memory kernel for the Gaussian chain model is calculated explicitly and a closed 
equation for the non-ergodicity parameter is derived. Finally, section V discusses the 
main results and perspectives. 

II. THE GENERATING FUNCTIONAL OF A POLYMER MELT IN TERMS 

OF THE COLLECTIVE VARIABLES 

Consider a homopolymer melt of M chains with the p— th chain configuration at a 
time moment t characterized by the vector function R( p )(s,t), where s numerates the 
segments of the chain, < s < N. The simultaneous dynamical evolution of R/ p )(s,t) is 
described by the Langevin equations 



5H{Rf} 




(1) 



with the Hamiltonian 






and the random force correlator 




(3) 



5 



where £o denotes the bare friction coefficient and from now on the Boltzmann constant 
k h = 1. 



After using the standard MSR-functional integral representation [21, 22 1 for the system 



TO, the GF takes the form 



3 M 



Z{lf\ s ,t),lf(s,t)} = / J] U DR<f\s,t)DRf ( a ,t)exp \a [R^ , £j 

j=ip=i { 

+ -YY dt dsds' -^ikdRf(s,t)V(k) exp \tk[R ip) {s,t) - K im) {s',t 
2 p=lm=l J Jo J ( 27r ) 

+ ]T / dt / cfe (s, t)/j p) (s, t) - iRf(s, t)i { f\s, t) 



(4) 



where the MSR-action of the free chain system is given by 



m r r N ( 
A {Rf\ Rf} = £ / dt / d J TG,[*R W (*. *)] 



<9 



5 2 



eo^( S ,t)-e^( S) t) 



(5) 



where £ = M- is the bare elastic modulus of a spring with the length of a Kuhn segment 



denoted by /, V(k) is the interaction energy of the chain segments, lj and 1^ are external 
fields conjugated to iffl and R^f 1 respectively and the summation over repeated Cartesian 



indices is implied. Now by the same way as in [23] collective variables can be introduced. 



As opposed to the statics [ITBHT^] we need to consider not only the mass density 



M rN 



p(r,t) = W dsi(r-RW(M)) 

i Jo 
p=i 



(6) 



but also the longitudinal projection of the response fields 



M rN 



U(r,t) = E / ^i J Rf ) (s,t)V i 5(r-R (p) (s,t)) 
p=i ® 



(7) 



where again the summation over repeated Cartesian indices is implied. Then the GF 



becomes 



. M 

Z{ - }= / l[DK^{s,t)DR^{s,t)DpDU 

P =i 



x 5 



M 



x 5 



p(r,t)-£ / ds5(r-Bp>\s,t)) 

M . 

n(r,t) - E / daikf(s,t)Vj5(T - R {p \s,t)) 
P =i J 

x exp|-i|^V 1 rfV 2 n(r 1 ,t)p(r 2 ,t)V r (r 1 -r 2 )+A {R (p) ,R (p) } 

where the dots imply some source fields which will be specified later. 
It is convenient to introduce the 2-dimensional field 



(9) 



where a = 0, 1 and 1 = (rx,£i) is used for abbreviation. In terms of the 2-dim. density 
(0) the GF (H) takes an especially compact form 



Z {V> Q } = / Dp tt (l) exp - - Pa (l)f/ Q/3 (12) P/3 (2) + W{p a } + p a (l)^«(l) (10) 



with the action of the free system 

. A/ 

W^{p,n} = ln / n^R (p) (s,t)^R. (p) (s,t)exp{Ao{R (p) ,R (p) }} 
J P =i 

M . 

p(v,t)-J2 / <Mr-RM(a,t)) 

Af . 



x 5 



and the 2x2-interaction matrix 



U af3 (l,2 



V(|n-r 2 ^ 
.n|ri-r 2 |) 



(12) 



and the ip a (l) is a source field conjugated to the 2-density @. In eq.([H]) and below 
the summation over repeated Greek indices and integration over variables with bars is 
implied. 



The exact form of the action V^jpo,} is not known explicitly, but can be obtained by a 
functional expansion by assuming that the density fluctuations are not very large and the 
functional iy{p a } is convex. 



The calculation is quite similar to the static case [lT6Hl^|. Let us introduce the cumulant 
GF, the "free energy", of the free system 



F{ijj a } = lnZ {tf; a } 

= ln J D Pa exp {W{p a } + Pa (l)M±)} 



This GF has the expansion 



FWa} = FfHlWaQ) + ^^(12)^(1)^(2) + 

+ ^^(123)^(1)^(2)^(3) + • . . 



with the free system cumulant correlators 



fP(i) 



^2(1,2) 



2, 3) 
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5 3 



1/1=0 



^(1)^(2)^(3) 



ip=0 



V>=o 



We are searching for the action H^{p Q } in the form of an expansion 



W{ Pa } = W{( Pa ) } + -W$(l2)5 Pa (l)5pp(2) + 



2! 



+ ^wi%(m)8p a (l)8pp(2)8 Pl (3) + 



3! 



where 



(13) 



(14) 



(15) 
(16) 
(17) 



8p a (l) = p a (l) - (p a (l)) 
= Pa(l)-F«(l) 



(19) 



In order to determine the coefficients in the expansion (|I8|) one should use the saddle 



point method when calculating the functional integral (fig). This can be carried out in 



the same spirit as in refs. [17,18|. This results in a Legendre transformation with respect 



to the extremum field p a {l)'- 



F{i; a } = W{p a } + p a (l)MK 



(20) 



As a result 241 it is found: 



P«(l) 



5F{ip a } 
SW{p a } 

Sp a (l) 



By the use of the expansions (|I4] ) and (|18|) in eqs. (^T 
with the free system cumulant correlators [H]Jz 



^(1,2) 



-1 



'1,21 



a/3 



^Ul,2,3) = F^(l,2,3) 



W^(1,2,3,4) = F^(1,2,3,4) 



(21) 
(22) 

we can relate the coefficients 



- F. 



(3) amp. 



af3~/ 



l,2,3)F^(33)FPi amp -(§,3,4) 



(3) amp. 
77<5 



(23) 
(24) 



(25) 



where the so called amputated correlators are defined by the relation [24]: 



w (n)amp. ( -, 9 \ 
*a/3-5 11)2, ...,nj 



) amp. 



l...,n) 



The main result of this section is the GF ([!(]) given now by 



Z{ip a } = / Dp a exp 



a/3 



;i, 2)^(1)^(2) 



+ ^^(123)^(1)^(2)^(3) 



+ ^^(1234)^(1)^(2)^(3)^(4) + 



(26) 



+ p a (l)^ a (l) 



(27) 



where the coefficients in the effective action are expressed in terms of free polymer system 



dynamics and given by eqs. (|23| - |25| , p6|) . The relation ( f27|) is the dynamical generalization 
of the coarse grained partition function which was obtained (for a diblock copolymer melt) 



in the ref. [17, IS 



III. THE EQUATIONS OF MOTION FOR THE TIME CORRELATION AND 

RESPONSE FUNCTION 

The representation of the GF (|27|) is a good starting point for taking into account 
fluctuation effects which enable us to go beyond the standard RPA methods. Before we 
proceed in this direction let us rederive the simple RPA-results for convenience and as a 
consistency check. 



A. RPA-results 



If we restrict ourselves in the expansion (|27|) to the quadratic order, then the RPA- 
correlator and response function S a p(l, 2) is obtained as a 2 x 2-matrix form 

S^(l 1 2)=[u + F^- i r 1 (l,2) (28) 

Here U is the matrix of interactions and given by eq. (|12|) . The correlation 2 x 2-matrix 
of the free system has the form 

F P) ri 2 x ( F oo(h 2)^(1,2) \ 

F ^ 2) ~ [fS\i,2) J (29) 

In eq. (j2~9"|) off diagonal elements FqI\i, 2), Fj$ (1, 2) are retarded and advanced suscepti- 
bilities respectively; the diagonal element F q (1, 2) is the density correlator. The relation 
between them is determined by the FDT in (k, ^-representation 
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Pj/St (k, = F$ (k, t) - F# (k, t) 



or alternatively in (k, ^-representation 

$iujF$(k,uj) = FS\k,u) - F (2) (k, W ) 



(30) 



(31) 



where f3 = 1/T. The inversion of the 2 x 2-matrix in eq. fl28|) results in the matrix elements 
in the (k, ^-representation 



S 00 {k,uj) 

S 01 (k,u) 
S w (k,u) 



1 + V(k)FS\k,u)} [l + V(k)F$\k,u>) 
Fg\k,uj) 



l + V(k)F£\k,u) 
l + V(k)F$(k,u) 



(32) 

(33) 
(34) 



From eqs. (|32| - |34]) we can simply see that the FDT for the non- interacting system ( p0| , |31|) 
assures the validity of FDT in RPA 



(3iuS 00 (k, lo) = S w (k, u) - S i(k, u) 



(35) 



Eqs. (|32| - |3"4"D coincide with the classical expressions for RPA-susceptibilities which have 
been used in the theory of liquids (see e.g. eq. (2.7.42) in |TT|| ) and more frequently in the 
dynamic theory of polymer melts P^-|T5|. 



B. Equations of motion beyond RPA: Role of fluctuations 

In this general case the full renormalized correlator 2 x 2— matrix G a p{l, 2) must satisfy 
the Dyson equation 





Jl,2) = 


s- 1 




ap 





a/3 



(1,2)-£ Q/3 (1,2) 



(36) 
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where the so called self-energy functional E Q!| g(l,2) contains only one-line irreducible di- 
agrams (or diagrams which cannot be disconnected by cutting only one line) |[24|| . These 
diagrams are built up from the vertices given by the effective action (|T0|) and lines with 
the full correlator matrix G Qi a(l,2) assigned to them. The diagrams which are relevant 
for the glass transition problem will be discussed later. 
The self-energy matrix has the triangular block structure 



S^(l,2) 







E„i(l,2) 



E 10 (l,2)E n (l,2) 



Inversion of the Dyson equation (^6|) yields the result 



G 00 (k,uj) 

G 01 (k,u) 
G w (k,uj) 



S o{k,u)S i(k,u)S 10 l (k,uj) + Su(k,w) 



S i(k,u) - Eio(k,w)J [S^^u) - E i(k,u;) 
1 



Soi (k,w) - Ei (k,cj) 



S^ik,^) - E i(k,o;) 



(37) 



(38) 

(39) 
(40) 



If besides relation ( 35|) FDT is also valid for the full correlator and response functions 



(3icvG 00 (k,uj) = G w (k,w) - G i(k,o;) 



(41) 



or in time domain 



(3^G 00 (k,t) = G 01 (k,t)-G 10 (k,t) 



(42) 



then the exact self-energy obeys 



PiuE n (k,u) = E i(k,o;) - Ei (k,a;) 



(43) 



or correspondingly in the time domain 



/3^E 11 (k,t) = E 10 (k,t)-E i(k,t) 



(44) 
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This agrees with the corresponding relations for self-energy given in [ 25] . We stress that 
the FDT for the free system (|3l|) as well as in RPA fl35[) is always valid. For the full 
correlator and response functions of a glass forming system this is not obligatory the case 



In order to proceed further let us note that for time intervals t > tr (where is the 
characteristic Rouse time) one can use for the free system correlator the diffusional ap- 
proximation pu 



F£\k,t) = F st (k)e 



-k 2 Dt 



(45) 



where D is the Rouse diffusion coefficient and F st {k) is the corresponding static correlator. 
Then taking into account FDT fl3T|) we have 



F$(k,u) 



Fg\k,u) 
i4 2) (k,c) 



2k 2 DF st {k) 
oo 2 + (k 2 D) 2 
f3k 2 DF st {k) 
-iu + k 2 D 
Pk 2 DF st (k) 
iuj + k 2 D 



(46) 
(47) 
(48) 



One should use eqs. ( f46"l - fi8| ) in the RPA-result (p2| - ^41 ), then the expressions for the full 
correlator and response functions (|38| - |40|) , after going back to the time domain, yields 



T 4 + *» (k) 

+ X7t 1 (k) 



(49) 



G i(k;M')+ / ^(k^^^Go^k;*",*')^"^-^-*') 
G 00 (k; t, t') + f E 10 (k; t, t")G 00 (k; t" , t')dt" 

J — oo 

+ f S 11 (k;t,t")G 10 (k;t",t')^" = -2Tr c G 10 (t,t') (50) 

J — oo 



where the inverse RPA-static susceptibility 



X ; t 1 (k) = [(3F st (k)]- 1 -V(k) 



(51) 



and the bare correlation time 
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Pk 2 DF st (k) 



(52) 



are defined. The initial conditions for the equations and ( EUJ ) has the form 



r c G i(k;t = t' + + ) = -1 



G i(k;t = O =0 



(53) 



and 



Goo(k;t = t / )=G', i (k) 



(54) 



where G st (k) is the full static correlator. Eqs. (^,^) represent the general result of the 
present paper, and below we will discuss the physical aspects of them in more detail at 
specific examples. 



The resulting equations (|49| - p4] ) are indeed very general. We made only use of the 
diffusional approximation for the free system correlator (H5|) and of a causality condition. 



Qualitatively the same equations of motion was obtained in the dynamical Hartree ap- 
proximation for a test chain in a melt (see eqs. (29)- (32) in fl2"0| ) and for a manifold in a 
random medium (see eqs. (B3)-(B5) in [0). Let us now consider the case when besides 
eqs. ( [4~9| - |5"4] ) the FDT (f42"D for the full correlator and response function holds. 



C. Time-homogeneity and FDT for the full matrix G a p{\L;t) hold 



Let us assume in eq. fl50|) t' = and t > 0. Then after differentiation of both sides of 
eq. (|50D with respect to the time and taking into account FDT (|42] ) we have 



TcJ^ + XrtOO G i(k;t) + j^ 10 {\i-t-t')G ol {\i-t')dt' 
+ dt' |/3 J^n(k; t-t')- E 10 (k; * - t') | G 10 (k; t') = 



(55) 
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The comparison of the eq. (|55f) with the eq. ( fl9|) yields 



d 



/3-S n (k;t) = S 10 (k;t) 



(56) 



which is again the familiar relation ([H]) for t > 0. As it should be the case, one of the 



eqs. (||], 0) is getting redundant now. 

Using the eqs. (f42[) and (HI) in eq. ([50]) and after integration by part we arrive at the 



result 



G o(k; t)+pf E u (k; * - t')-^G 00 (k; t')dt' = 



(57) 



where 



x; t 1 (k)= [3F$\k) --y(k)-/3E u (k ; * = o) 



(5f 



The last term in eq. (|58f) is the contribution of fluctuations in the static correlation 
function. It must be stressed that eq. fl57D is equivalent to the Mori-Zwanzig equation, 
derived by the projection formalism || and was specified for the two slow variables, 
mass density and longitudinal momentum density, in ref. |L|. Performing the Laplace- 



transformation 



£(•••) = / dt . . . exp(izt) 
Jo 



(59) 



in eq. ([57]) for 4>(k, z) = G 00 (k, z)/G st (k) we get 



(k,z) 



1 



. fc 2 G;/(k) 

r + fc 2 M(M) 



(60) 



where 



To 



DF st (k) 
M(k,z)=(3 2 E u (k,z) 



(61) 
(62) 
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and we have used the relation 

Xjt 1 = TG-\k) (63) 

The interesting point is that the self-energy matrix element En is proportional to the 
memory-kernel. Indeed eq. ( |62| ) connects the Mori-Zwanzig technique with the MSR- 
formalism. The particular form of the self-energy matrix element En for the physical 
problem defined by the action given by eq. (|I0| ) will be discussed in the section IV. 

It can be expected that a critical temperature T c exists where the correlator Goo(k, z) 
as well as the memory kernel acquire a pole at z = 0. This would show an ergodicity 
breaking or an ideal glass transition [0 and will be considered in section IV. 



D. The time- homogeneity is valid but FDT is violated 



This case was discussed in the literature [26|,30,31]. The more general case, when the 



time-homogeneity does not hold any more was also considered [27,28 



According the scenario given in [p6| , pO| , pT| below the temperature T c of the ergodicity 
breaking the phase space decomposes into regions (ergodic components). The latter are 
separated by high barriers which can not be crossed at times t < t*. The time t* has 
the physical meaning of a "trapping" time. Especially for a system with infinite range 
interactions, we have t* — > oo in the thermodynamic limit. This point had become most 
obvious in the case of spin glasses with long range interactions [p6|l . In such systems 
quenched disorder is present from the beginning and it has been shown that replica sym- 
metry breaking corresponds to the splitting of the phase space by infinite high barriers. 
The present dynamical theory on conventional glasses without quenched disorder sug- 
gests a similar scenario. Disorder develops during cooling and eventually large barriers 
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develop. Then these barriers cannot be crossed and the phase space is broken up into 
accessible parts. Dynamics becomes slow and glassy. Mathematically this is monitored 
in the invalidity of the FDT. 

To be more precise, we assume that for our case t* is large but still finite. For the 
intervals t << t* the relaxation occurs within one ergodic component, the dynamics is 
ergodic and FDT holds. For t » t* the system jumps over the phase space barriers 
from one component to another and FDT is violated. In this case the total response 
to an external field consists of two parts. One of them is determined by the dynamics 
inside one ergodic component (intracomponent dynamics) obeying FDT, and the other 
which appears at t > t* as a result of crossing the barriers between different components 
(inter-component dynamics). The latter process violates the usual FDT. 
With this in mind let us make the following assumption 
d 



ft 



^+ T sign(t)0(|t|-|**} 



Goo(k,t) = G 01 (k,t)-Gi (k,t) (64) 



which might be the linear version of a more general relation 



^ + 7 X{G 00 (k,t)} 



Goo(k,t) = Goi(k,t)-Gi (k,t) (65) 



In eqs. (p4] , |65] ) the (phenomenological) parameter 7 has the dimension of an inverse time 
and, as we will see below, has the meaning of a characteristic rate of the inter-component 
dynamics ("hopping process" in the nomenclature of Ref. ||). X{G o} is an arbitrary 
functional of the correlator G 00 (k, t). In eq. (|64"D #(•••) is the ^-function and sign(t) 
keeps the correct transformation under time reversal: Goi( — t) = G\o(t). We will call 



eqs. ©H) after ref. ggfflj Quasi-FDT (QFDT). 

Let us use eq. |>§) (at t > t* > 0) in eq. (0) ( at f = and t > 0) by acting with 
(3(8/ dt + 7) on its both sides. After the same calculations carried out already in section 
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Ill C we arrive at 



P^+l) £ii(k;*) = £ 10 (k;f) 



(66) 

where t > t* > 0. As before, the eq. ( |66"D assures that one of the eqs. gets 
redundant, whenever the QFDT flB^D is valid. 

The substitution of eqs. ( |6lD and (|66|) into eq. (|50|) yields 



d_ 



.—1 



T c^ + Xst 



Coo (k; + £ E u (k; t - ^G 00 (k; dt 



+2/5 7 / * Eu(k; t - t')G 00 (k; t')dt' + £7 /' * S u (k; t - i')Goo(k; = (67) 
Two important limiting cases can be distinguished: 

• If t — *■ 00 and £* — ► 00, but t/t* — > 0, then the two last terms in eq. ( |57|) can be 
neglected and we go back to eq. (|50|), which contains (under conditions discussed 
in section IV) an ideal glass transition. The system could never escape from one 
ergodic component (an absolutely confined component in nomenclature of ref. ||32|| ). 

• If t — > 00 and t* is large but finite, (so that t*/t — > 0), then the next to last term in 
the l.h.s. of eq. ( p7| ) can be neglected. In this case the eq. ( 57 ) takes the form 



Tc dt +Xst 



G 00 (k; t)+/3 J* E u (k; t - t') (JL + 7) G 00 (k; t')dt' 



(6* 



The Laplace transformation of eq. flS8|) gives the result 



<i>(k,z) 



-iz + 1 + - 



k 2 Xst(k)/T 



(69) 



T° + k 2 M(k, Z ) 



with 



Xst = Xs^ - l T c 



(70) 
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In contrast to the ideal glass transition case at z — > the kernel M(k, z) is large 
but finite. Then , according eq. ( |69"D at z — ► the behavior of 0(k, z) is determined 
mainly by a simple pole at = 7. Instead of going to a plateau the correlator 
0(k, i) decays with the characteristic rate 7. This is a result of the inter-component 
dynamics. 



IV. THE EXPLICIT FORM FOR THE MEMORY KERNEL M(K,Z). THE 
MODE-COUPLING APPROXIMATION. 

For the ideal glass transition problem fl],§] the correlator Goo(k, z) as well as the 
memory kernel M(k, z) acquire a pole at z = at the critical temperature T c . One 
can easily see that for the effective action, given by eq. (|T0|), such contribution come 
from the sum of all "water-melon" diagrams which are represented in Fig.l. Each line 
denotes the full matrix G a/ 3(k, z) and a vertex with m legs denotes the bare vertex function 
wffi ~(ki, Z\\ k 2 , z 2 ; ■ ■ ■ k m , 2 m ). In the MCA all vertices renormalization is neglected |33| . 



Another type of diagram, the "tadpole" diagrams, which are shown in Fig.2, appears in 
the context of a Hartree-approximation |30||3T||34]1 . On the other hand the contributions 



of these diagrams remain finite at z — > and because of this they are not relevant 
for the ideal glass transition problem. On the contrary, these diagrams are essential in 
the context of e.g. fluctuation effects in the theory of microphase separation in block 
copolymers [|l8l , |l9f . It was shown there, that fluctuations change the order of the phase 



transition from two (mean field) to one (upon renormalization) . In this research field the 
corresponding procedure has been terminated as Brazovskii renormalization. 

The explicit expression for the arbitrary vertex function W^ Jl, 2, ... n) is not known 
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in detail. That is why we can restrict ourselves only to the first diagram given in Fig.l, 
which corresponds to the one-loop approximation. As a result we have 

Sn(k; z) = 2(i) 2 J ^w[%{k, k - q, q)G m (k, k - q)G >y(g)W§> (g, k — q,k) (71) 

where the short hand notation, k = (k, z) and q = (q, s), was used and the expression for 
the vertex function has the specific form 

1^(1,2,3) = F^(l,2,3)Fg- 1 (l > l)l^- 1 (2,2)F« ) - 1 (3 > 3) (72) 

In general there are nine terms in the sum (|7T|), but only one of them, which does not 
include a response function has the important l/z-singularity at z — > 0. This leads to the 
expression 

S n (k, z) = 2(^) 2 | -0^w£l(k, k - q, q)G 00 (k - g)G 00 (g)^ ( 3 J (q, k - q,k) (73) 
where 

Wjj&ik, k-q,q) = F®(k, k - q, q)(F ( - 2 ' ) ~ 1 ) 01 (k)(F ( - 2) ' l ) w (k - ^(i^- 1 )^) (74) 
and 

W^(q, k-q,k) = F®(q, k - q, k)(F ( - 2 ' ) ~ 1 ) 1 o(q)(F ( - 2 ' ) ~ 1 ) 10 (k - q^F^Mk) (75) 

In eq. (ff3|) the integrals over Laplace-frequency s is taken along the straight line in the 
complex s-plane above all singularities of the integrand. Since a pole of Goo(q, s) at s=0 
predetermines the l/z-behavior of the whole integral, we are able to consider the vertex 

(3) (3) 

functions and Wool om y i n the static limit, s — > and 2 — > 0. Below, we give these 
limits for the 2- and 3-point response functions. 

The Laplace transformation of the 3-point response function, which appears in eq. flT3|), 
is given by eq. (|A(]) in the Appendix. As a result its static limit reads as 
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lim Fg\(k 2 z 2] k 3 , z 3 ) = 2 F£>(k 2 , k 3 ) 

22,23^0 



?2c,(3). 



(76) 



The static limit of the 3-point response function, which appears in eq. (|75|) has the same 
form 



lim F£l(z 1 ,z 2 )=FF£\k 1 ,k 2 ) 

22,23^0 



?2c(3). 



(77) 



The static limit for the 2-point response functions is given by the relation 



limF^k,*) = limF^(k^) = -(3F st (k) 

z—>[) z— >0 



^(2), 



Explicitly the expressions for F$> and are given by [[36 



(7S 



fS\k u k 2 ) 



-N 2 p 



with 



J a (K?) - J 2 (Ki) , J 2 (Ki) - J a (K§) J a (K§)-J 2 (K?; 



K?-Kg 



+ 



K|-KI 



K§-K? 



(79) 



Ki + K 2 + K 3 = 



(80) 



and 



f£\k) = n Pq j 2 {K' 



(81) 



where K = kL/N/Q and p = c ^ is the average segments concentration. The function 



Ux) = 2 



e x — 1 + x 



(82) 



is known as the De-bye function in polymer physics and is usually approximated well by 
the more simple Pade expression P^,f3?J 



Jo.(x) 



1 + \x 



(83) 
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By making use of eqs. (|79| - |83|) the expression for the vertices (ffi|, ff5D in the static limit 
becomes 



iyi t 3) (K 1; K 2 ,K s 



1 



A(3piN 



3 + -(K? + K2 + Kl) 



34) 



where K 2 = k 2 Z 2 iV/6. Then for the memory kernel (62) we derive 

d 3 qds 



.Ufk.:)=2(if-i 



A\J piN 2 J (2vr) 4 
x k — q, — z — s)0(q, s 



3 + ^V + g 2 + k q ) 
o 



G st (-k-q)G st (q) 



15) 



The integral over q is mainly determined by the strong peak of G st (q) at q = qo = 1/er, 
where cr is the bead diameter in the spring-bead model for the chains in the melt. In 
that case the second term in the brackets dominates and iV-dependence for long chains is 
cancelled, as it should be. 

In the glass state the correlator 0(k, z) has the form 



lim#k,z) = ^ 



(86) 



where /(k) is the non-ergodicity parameter. From eqs. |6^ and (p5|) one can easily see that 
the function /(k) satisfies the equation 



^[k 2 + q 2 + kq] 2 G st (-k - q)G st (q)/(-k - q)/(q) 



57) 



which qualitatively corresponds to the result of the conventional MCT (see e.g. eq. (3.37) 
in ref. [|l[). The factor Z 4 in front of the mode-coupling integral (j87D indicates that the 



relevant length scale for the glass area is indeed the Kuhn segment length. This confirms 
statements, that in polymer melts, where a large variety and wide range of internal degrees 
of freedom dominate the physical behavior the glass transition is indeed ruled on local 



22 



scales, i.e., the range of the nearest neighbor interactions. Moreover, the present model 
suggests, that the glass area is larger for stiffer chains, since I becomes larger. This is 
qualitatively correct but needs a more detailed investigation. 



V. SUMMARY AND DISCUSSION 

In this paper we have shown that the MSR-functional integral representation is very 
convenient for the treatment of the Langevin dynamics of the polymer melt in terms of 
the collective variables: the mass density and the response field density. The expansion 
of the free energy and the Legendre transformation technique, which was given in the 



ref. []16| |18[1 for the static case, was extended here for the dynamics. As a result we have 



derived the dynamical GF (0) with the action that allow results beyond the RPA up to 
arbitrary order of the density and/or response fields density. 

It was shown that the GF is a good starting point for the derivation of the general 
dynamical equations. This is the set of equations ( fi5| , pD| ) for the correlation and response 
functions, with the free part determined by RPA. For the particular case when the time- 



homogeneity and FDT are satisfied these two equations reduce to the one eq. (|57|). 



It is obvious that eq. ( [57|) is equivalent to the Mori-Zwanzig equation with the memory 
kernel given by the matrix element Sn(k, t). In the framework of the Mori-Zwanzig 
formalism and MCA the glass transition problem was extensively discussed in In 
our approach these results are qualitatively restored if the "water-melon" diagrams for 
Sxi(k, z) (see Fig.l) are summarized. The results of the vertices renormalization could 



be in principle also investigated f3"Tfl. 



The situation when FDT is violated is much more obscure in spite of the extensive 
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discussion [^6H28U30|j3l|| . It could be guessed that the physical reason for this (at least 
at a finite £*) is an intercomponent dynamics as it was discussed in section IV D. We 
have made a simple but plausible assumption about the form of QFDT (0). When 
this assumption was made in the corresponding model calculation we had shown that a 
sharp (or ideal) glass transition is smeared out by this intercomponent dynamics with a 
characteristic rate 7 and the density correlator takes the form (|69|). 

For the ideal glass transition case we have calculated the expression for the memory 
kernel (BH) by using the explicit forms for the free system's static correlators. As a result 



the equation (|87]) for the non-ergodicity parameter f(k) is qualitatively the same like 
in MCA. The physical picture we had derived here is most interesting. We have indeed 
shown that the violation of the fluctuation dissipation theorem yields a similar picture as 
in systems where quenched disorder is present from the beginning. For example in spin 
glasses it was shown, that the presence of disorder and frustration yields a glassy phase. 
Moreover it could be shown that the case of replica symmetry breaking is responsible 
for the disordered nature of the phase space, especially for the large barriers dividing 
the phase space. Thus the entire phases space is no longer accessible for the system. 
This corresponds e.g. to glassy dynamics. In structural glasses such as polymer melts 
which form easily glasses, quenched disorder is not present. In the melt phase equilibrium 
dynamics determines the structural properties, even though the dynamics is very slow 
mainly ruled by the molecular weight. Here we must start from equations that take into 
account all interactions. With the formalism presented in the sections above we succeeded 
to derive a similar picture on totally different physical grounds. We had shown that 
glassy phases could appear upon violation of the FDT. At a simple model computation 
we had presented arguments, that a characteristic time t* can be related to barrier heights 
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that eventually can be crossed. This corresponds to the breaking of the phase space into 
different regions that are separated by barriers, which have their origin in the interactions. 
In fact, these barriers correspond to the development of a multi valley structure of the 
phase space. We will come back to this point in later publications. 

The present general formalism allows important generalizations for interacting sys- 
tems. For example we will treat homopolymer blends and copolymer melts which tremen- 
dously enriches the picture of the dynamical behavior. The dynamics of copolymers was 
studied above a microphase separation temperature T ms in RPA [PLIKJ and below T ms 



by the numerical solution of the Ginzburg-Landau equation [[38 j. One can easily derive 
a copolymer counterpart of the dynamical GF fllCf ). Then the Hartree approximation, or 
summation of all tadpole diagrams shown in Fig. 2 (in the same manner as for the static 



18| . |T9|| ) gives a direct way to obtain a closed dynamical equation for the composition- 



composition correlation function. Moreover we have now the possibility to study the 
glass transition in blends and copolymer melts. These problems are of wide experimental 
interest. Imagine for example that one component of the blend or one species of the 
copolymer melt freezes out during cooling. The striking problem is then to discuss the 
interplay between freezing and phase - or microphase separation. The dynamics of one 
species becomes very slow and some parts of the systems become eventually immobile 
at certain correlation length, corresponding to the distance of the (ideal) critical point. 
These problems are also under current investigation. 
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APPENDIX A: THE 3-POINT RESPONSE FUNCTION AND ITS LAPLACE 

TRANSFORMATION 



The 3-point response function which appears in eq. ([71) is determined by 
F^l(k 1 ,ti,k 2 ,t2-,k 3 ,t 3 ) = (ik 2 )j{ik 3 )i / ds 1 ds 2 ds 3 (Al) 

Pl,P2,P3 = l ° 

Rf 2) (s 2 , t 2 )iRl m) (s 3 , t 3 ) exp {ikiR (pi) (si , t x ) + ik 2 R {p2) (s 2 , t 2 ) + ik 3 R> 3) (s 3 , t 3 )}) Q 



x (i 

where < • • • >o stands for the averaging with the action of the free chain system. By 
making use of the Nonlinear-FDT (NFDT) rule (see eq. (2.31) in ref. ||35|| ) 

-iRj(s,t) -^p^-Rj( s ,t) (A2) 

we will come to the following relation 

Q 2 

F ?l(ki, *i; k 2 , t 2 - k 3 , t 3 ) = /3 2 ^-FS(ki, ti; k 2 , t 2 ; k 3 , t 3 ) (A3) 
with the causality condition t\ > {t 2 ,t 3 }. In the same way we get 

Q 2 

Fg(k 1; t r , k 2 , t 2 ; k 3 , t 3 ) = [3 2 -^-F$ Q (k 1; t x] k 2 , t 2 ; k 3 , t 3 ) (A4) 
with the causality condition t 3 > {ti,t 2 }. 

Let us define the Laplace-transformation of the 3-point response function 
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(3) f°° f f (3) 

Fqu(zi, z 2 , z 3 ) = / dtt / dt 2 dt 3 F$ 1 [(ti,t2,t 3 )exp{iz 1 ti+iz 2 t 2 + iz 3 t 3 } 

J —oc J —oo J— CO 

/•o /-o (3) 

= 2tt5(zi + z 2 + 23) / dt 21 / d* 3 i^oii(*2i, *3i) exp{z«2*ai + i*3*3i} 

= 2tt5(*i + z 2 + z 3 )f£\(z 2 , z 3 ) (A5) 

where we have used the causality condition t\ > {t 2 , t 3 } and the time translational invari- 
ance. 

By making use of the NFDT ( |A4| ) and after integrations by parts we have 

F ( g(k 2 , z 2 ; kg, z 3 ) = (3 2 ^FS\k 2 , k 3 ) - ^ 2 F ( 3 J(k 2 , z 2 ; k 3 , t 31 = 0) 

- ^3^ooo( k 2, hi = 0; k 3 , z 3 ) + {iz 2 ){iz 3 )F^{k 2 , z 2 ; k 3 , 2 3 ) j (A6) 

where 

F ( oJ(k 2 , z 2 ; k 3 , z 3 ) = ( dt 31 [ dt 21 F^l(k 2 , t 21 ; k 3 , £ 31 ) exp{^ 2 t 21 + iz 3 t 31 } (A7) 
J —00 J— 00 



and 



F ( oJ(k 2 , 2 2 ; k 3 , £ 31 = 0) = f dt 21 F$ (k 2 , t 21 ; k 3 , i 31 = 0) exp{?,2 2 t 2 i} (At 

•/ —00 



and Fst\~k 2} k 3 ) is the static 3-point density correlator. From eq. (|A6|) we immediately 
obtain the static limit (|76| ). 
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FIGURES 




?IG. 1. Diagrammatic representation of the self-energy En (k, z) in MCA, which has a simple 
pole at z = 0, i.e. is relevant for the ideal glass transition. The vertices are bare: W$\W^\..., 
i.e. this approximation neglects all vertex renormalization. 



31 



W 



7IG. 2. The sum of all tadpole diagrams, which are finite at the external frequency z — ► 0. 
They are not relevant for the glass transition problem. 
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